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Charge fluctuations for particles on a surface exposed to plasma 
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Abstract 

We develop a stochastic model for the charge fluctuations on a microscopic dust particle rest- 
ing on a surface exposed to plasma. We find in steady state that the fluctuations are normally 
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distributed with a standard deviation that is proportional to (CT e ) ' , where C is the particle- 
surface capacitance and T e is the plasma electron temperature. The time for an initially uncharged 
ensemble of particles to reach the steady state distribution is directly proportional to CT e . 
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Figure 1: Schematic of plasma assisted electrostatic cleaning (PAEC). A particle on a surface 
exposed to plasma acquires a charge Q < and is pulled from the surface when the electrostatic 
force QE W due to the sheath electric field E w exceeds the adhesive van der Waals force -F v aw- 

As features on integrated circuits become smaller, the problem of contamination by mi- 
croscopic particles becomes larger. The best solution is to minimize particle creation and 
to prevent particles from reaching the substrate. However, if particles are deposited on 
the substrate then methods for removing them are needed. One proposed method is called 
plasma-assisted electrostatic cleaning (PAEC). 

PAEC works by exposing a surface with dust particles on it to a plasma (Fig. [TJ. It is 
conjectured [I] that a particle can acquire a net charge Q < from the fluxes of charged 
plasma species. The electrostatic force on the particle due to the sheath electric field E w 
is QE W . When this force exceeds the adhesive force the particle is pulled from the surface 
and experiences a large acceleration across the sheath. If there are no nearby regions where 
particles can be trapped [21 E], the particle will move away, and the surface will have been 
"cleaned." Thus, the plasma serves two purposes: charging the contamination particle 
and enhancing the electric field at the substrate. Several experiments [H IU [5] have shown 
that plasma exposure can remove micrometer to nanometer-diameter particles from both 
insulating and conducting surfaces. 

Sheridan and Goree [T] coated an aluminum surface with 0.2-10 /im alumina grains and 
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found that significant particle removal occurred in a low-density plasma (n e ~ 10 8 cm~ 3 ) 
when the surface was exposed to energetic electrons (59 eV) from an emissive filament. They 
observed that the removal rate increased with plasma density and decayed exponentially 
with time. Flanagan and Goree [I] used JSC-1 (Johnson Space Center lunar simulant one) 
particles on a glass substrate, again for low plasma densities and with energetic electrons, 
and verified that the time constant for removal decreases with increasing plasma density. 
They found time constants ranging from 30 s to Is for plasma densities from 2 x 10 6 cm -3 
to 2 x 10 8 cm -3 . Lytle et al. [S] performed experiments using a high density pulsed helicon 
source to remove particles as small as 30 nm from a dielectric substrate. Their results 
demonstrate that a minimum plasma exposure time is required for particle removal. In 
Lytle's studies the reported electron temperature was low (~ 3 eV), but the plasma density 
was 3 to 4 orders of magnitude larger than that in the other experiments [U H] . 

The dominant force holding a microscopic particle to a surface is thought to be the 
London- van der Waals force (HI |7j. The adhesive force on 10-100-nm diameter particles is 
reported [7] to be ~ 1-10 nN. An estimate for the van der Waals force between a spherical 
particle of radius a and a smooth, flat surface is [I] F v dw — Ha/(6h 2 ), where the Hamaker 
constant H ~ 10 -19 J and h is the separation between the particle and the surface (Fig. 
[I]). For a characteristic separation h = 0.3 nm and radius a = 10 nm, F v( ±w ~ 2 nN. 
Flanagan |4] estimated for his experiment that E w may have been as large as 10 kV/m . For 
this value of E w a particle charge Q ~ — 10 6 e would be required to overcome the predicted 
adhesive force. There is no evidence that microscopic particles can attain such large charges; 
the potential of a 10-nm radius particle in free space with this charge would be — 150 kV. 
Consequently, it seems likely that adhesive forces are significantly overestimated for typical 
processing surface conditions and that most particles removed by plasma cleaning are loosely 
bound. This may be due to surface contamination and/or to the fact that most particles 
are not perfect spheres and so rest on the surface at a small number of contact points whose 
asperity size is much less than the particle size, giving a reduced effective radius. 

One way that a particle on a biased surface can acquire charge is by sharing charge with 
that surface. Flanagan |4] referred to this as the "shared charge model." The shared charge 
on a spherical particle on a flat surface with electric field E w is [8] Q — (1.64)47reo-E' M ,a 2 , 
where the factor of 1.64 is a small enhancement because the particle is a perturbation to 
the otherwise flat surface. For microscopic particles under typical plasma conditions, the 
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shared charge is small. For example, for E w = lOkV/m and a = 500 nm, Q = — 2.9e. That 
is, on average, microscopic particles on a surface have \Q\ < e due to shared charge. Such 
an average charge is much too small to explain plasma assisted cleaning. 

Particles can also acquire charge because of electron and ion fluxes from the plasma. 
In steady state, these fluxes balance. However, fluctuations in the particle charge may be 
significant due to the discrete nature of the charging process. It is found experimentally [UH] 
that the rate of particle release during plasma cleaning decays exponentially, indicating that 
each particle release is an independent random event. This is consistent with release being 
due to particle charge fluctuations, where release occurs during large negative excursions in 
charge. In what follows, we develop a model of charge fluctuations caused by electron and 
ion fluxes to a particle on a surface. 

The ion flux across the sheath edge is Tj = uqc s (Fig. [T]) where uq is the plasma density 
at the sheath edge, c s = ^JeT e /mi is the ion acoustic speed and T e is the effective electron 
temperature in eV. We assume a single ion species with charge +e, and neglect secondary 
electron emission. Since the sheath is essentially source-free, continuity requires that Tj is 
also the ion flux at the substrate. The frequency with which ions strike a particle on the 
surface with ion collection area A{ is V{ = Y^Ai. The probability that a particle collects an 
ion in a time interval At <C u^ 1 is then Pi = z^At, which is assumed independent of the 
particle's potential. 

For electrons in thermal equilibrium (i.e., Boltzmann electrons), the electron flux to a 
surface at a potential with respect to the sheath edge is 

T e (0) = ^n ve^ T % (1) 

where v is the average electron speed. The average frequency with which electrons strike a 
particle with electron collecting area A e is then v e = T e A e , and the probability of collecting 
an electron in At is P e (0) = v e At. In the typical case where Tj < r e (0) the particle will have 
a steady-state floating potential 0/ < such that the electron and ion collection rates are, 
on average, equal. Assuming equal collection areas, A^ = A e , <f>f is a solution of z/j = v e (0/). 

Now consider a fluctuation around the floating potential, where the particle and surface 
are assumed to have the same average <f>f. The particle's potential can be written as <fi = 
<pf + 5<f), where 5<p is the potential difference between the particle and the surface. The 
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probability of electron collection then becomes 

P e = P^ T % (2) 

which decreases when 5<p < because electrons are repelled, and increases when 5(p > 
because electrons are attracted. This mechanism effectively limits the particle's charge 
fluctuations to a finite range about an average value S(f> ave = 0. 

Potential fluctuations can be related to charge fluctuations 5Q by 5(f) = 5Q/C, where C 
is the capacitance of the particle-surface system. The probability of electron collection is 
then 

Pe = Pi /Q/(CT e )_ (3) 

The charge on each particle performs a bounded random walk where the probability of 
gaining an ion (a step to the "right") is constant, and the probability of gaining an electron 
(a step to the "left") is given by Eq. The charging model [Eq. ([3f] represents a Markov 
process [10J for transition probabilities from a charge state 5Qj = je, j G (—00,00). Since 
the Markov process is regular and effectively finite the steady-state distribution of 5Q is a 
normal (i.e., Gaussian) distribution. This is our first result. 

Our second result is that the steady state, root-mean-squared value of the charge fluctu- 
ation (i.e., the standard deviation of the 5Q distribution) is 



SQrms / CT e 

— = y— (4) 

which gives the fluctuations in units of e. Consequently, SQ rms depends only on CT e . Note 
that SQ rms is not given by Poisson statistics for the net number of elementary charges on 
the dust particle [5], since in our model SQ^e <C SQ Tias . We also find that the steady-state 
value of SQ rms is independent of the plasma density, as is the case for isolated particles in 
plasma [TT]. To further characterize the charge fluctuations, we must consider C . 

The capacitance of an isolated spherical particle with radius a is Co = 47reoa, which is 
valid in plasma when a is small compared to the Debye length. If the particle is placed a 
height h < 0.1a above a flat conducting plate (Fig. [T]), the capacitance increases to [12] 

C 1 , 1, h 

_„ 7+5 h2--ln ; >l, (5) 

where 7 = 0.5772 ... is the Euler-Masherone constant. Here C diverges logarithmically 
as h/a — > 0, so that C > C , increasing the range of the charge fluctuations [Eq. Q]. 



The capacitance is insensitive to the shape of compact particles when a is taken to be the 
characteristic particle size [UJ. For values typical of Flanagan's experiment [1], a — 500 nm, 
h = 0.3 nm, and T e = 60 eV (the primary electron energy), we find C/Cq = 4.63 and 
<5Qrms = 3 lie, which greatly exceeds the shared charge. 

To determine the temporal behavior of SQ, the model was solved using a Monte Carlo 
simulation for an ensemble of n identical particles. For each particle, we start with the initial 
condition SQ = 0. For each time step At, we compute two random numbers distributed 
uniformly in [0,1). If the first random number is less than Pj, then SQ is increased by 
one, representing ion collection. If the second random number is less than P e , then SQ is 
decreased by one, represented electron collection. 

For a given value of CT e , we simulated this system for n = 10 000 particles vs dimension- 
less time tvi. During the course of the simulation the average charge SQ ave , the standard 
deviation SQ rms , the largest positive charge SQ max , and largest negative charge SQ min were 
computed for the ensemble. A time history of one simulation run is shown in Fig. [2] for 
CT e = 5.56 x 10~ 18 FeV, which corresponds to the capacitance of an isolated sphere with 
a = 20 nm for T e = 1 eV. Here the average of the charge fluctuations is approximately zero, 
while the extreme values make excursions of up to ±25e away from zero. When a particle's 
charge makes a large negative excursion, the electrostatic force may exceed the adhesive 
force leading to particle release. 

The dependence of the dimensionless charging time t r Vi on CT e computed using Monte 
Carlo simulations is shown in Fig. |3j Here t r is the characteristic time required for an 
initially uncharged ensemble to approach steady-state. The charging time was found by 
fitting SQ Tias (t) with a function oc (l — e~ t ^ tr ). Our third result is that the charging time 
increases linearly with CT e and is well described by the line 

Uvi = (2.05 x lO^F^eV" 1 ) CT e . (6) 

The charging time in seconds is 

t r = (72.6 x lO^eV- 1 ) oc — , (7) 

v an c s an 

where a is in m, no is in m -3 and c s is in m/s. Here t r is inversely proportional to the particle 
radius and the plasma density. The inverse dependence on plasma density is consistent with 
the measured jl] dependence of the particle release time constant. Smaller particles have a 
longer charging time since they have smaller collecting areas. 
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Figure 2: Time dependence of charge for an ensemble of n = 10 000 particles with CT e = 5.56 x 
10 _18 F eV. The dashed line gives the asymptotic dependence of Eq. &ty . 
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Figure 3: Dimensionless charging time t r Vi vs particle-surface capacitance C times the electron 
temperature T e in eV, determined from Monte Carlo simulations. The data are well fitted by a 
straight line. 
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Finally, we can roughly estimate the electrostatic force for a dust particle on an electrically 
floating surface in a discharge with a bulk component and a low-density energetic electron 
component [US]. For an electrically floating (or dielectric) substrate, the electric field at the 
surface scales as E w ~ T e ^ lot /\o, where T e j lot is the effective temperature of the energetic 
electrons and determines the floating potential. For a low-pressure filament discharge, T e ^ot 

1 /2 1/2 

is roughly the primary electron energy. The Debye shielding length ~ T e bulk n is 
dominated by the bulk electron space charge. Consequently, we predict the electrostatic 
cleaning force scales as 

F ex 5Q rms E w ~ a^XZtT^l (8) 

which gives a weak scaling with density and radius and a stronger scaling with the the 
tail electron energy. This scaling is supported by the experimental observations [H H] that 
particle removal occurs even at low plasma density when there is an energetic electron 
component. Since the adhesive van der Waals force oc a, it goes to zero faster than the the 
electrostatic removal force [Eq. ^8j)]. Consequently, there should be no lower limit on the 
size of particles that can be removed using plasma assisted cleaning. 
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